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Abstract 

The model under consideration is the two-dimensional (2D) one-com- 
ponent plasma of pointlike charged particles in a uniform neutralizing 
background, interacting through the logarithmic Coulomb interaction. 
Classical equilibrium statistical mechanics is studied by non-traditional 
means. The question of the potential integrability (exact solvability) of 
the plasma is investigated, first at arbitrary coupling constant F via an 
equivalent 2D Euclidean-field theory, and then at the specific values of 
r = 2*integer via an equivalent ID fermionic model. The answer to the 
question in the title is that there is strong evidence for the model be- 
ing not exactly solvable at arbitrary F but becoming exactly solvable at 
F = 2*integer. As a by-product of the developed formalism, the gauge 
invariance of the plasma is proven at the free-fermion point F = 2; the 
related mathematical peculiarity is the exact inversion of a class of infinite- 
dimensional matrices. 
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1 INTRODUCTION 



In this paper, we consider a classical (i.e. non-quantum) model which belongs to 
the general class of two-dimensional (2D) Coulomb systems of charged particles. 
According to the laws of 2D electrostatics, the particles can be thought of as 
infinitely long charged lines which are perpendicular to the confining surface. 
Thus, the electrostatic potential w at a point r, induced by a unit charge at the 
origin, is given by the 2D Poisson equation 

Av{r) = -27r5(r) (1.1) 

In a plane, the solution of this equation, subject to the boundary condition 
Vv{r) — > as |r| — > oo, reads 

^;(r) = -ln(M^ (1.2) 

where the free length constant ro, which fixes the zero point of the potential, 
will be set for simplicity to unity. In the Fourier space, the Coulomb potential 
(1.2) exhibits the characteristic small-A; behavior ?)(k) = l/|kp. This maintains 
many generic properties (like screening and the related sum rules [1]) of "real" 
3D Coulomb fluids with the interaction potential v{r) = l/|r|, r G R^. The pair 
interaction energy of particles with charges q and q' , localized at the respective 
positions r and r', is 

viv,q;T',q')^qq'v{\v-T'\) (1.3) 

A given continuous Coulomb system is classified via the number M of differ- 
ent mobile (pointlike) species a = 1, 2, ... , M, with the corresponding charges 
qa and particle densities Ua, embedded in a fixed uniform background of charge 
density p^. The most studied versions are the one-component plasma (OCP), or 
jellium, and the symmetric two-component plasma (TCP), sometimes called the 
Coulomb gas. In the OCP there is only one mobile species, M = 1, with q\= q 
and 111 = n, and neutralizing background of charge density pb. It is useful to 
introduce the background "number density" nb such that pb = —qrib; the neu- 
trality condition is then equivalent to n = rib- The symmetric TCP corresponds 
to M = 2, namely qi = q and ni = n/2, q2 = —q and n2 = n/2 (n stands 
for the total particles density), with no background charge pb = 0. Due to the 
logarithmic nature of the interaction, the equilibrium statistical mechanics of 
the underlying 2D Coulomb systems at the inverse temperature /? = l/{kBT) 
depends exclusively on the dimensionless coupling constant T = j3q^] the par- 
ticle density n only scales appropriately the distance. Both OCP and TCP are 
solvable in the high-temperature Debye-Hxickel limit F ^ (in the bulk [2] 
as well as for finite systems [3, 4]) and at the free-fermion point F = 2 (see 
reviews [5, 6]). Through a simple scaling argument, the density derivatives of 
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the Helmholtz free energy can all be calculated exactly at arbitrary F. For in- 
stance, the exact equation of state for the pressure P, /3P = n(l — r/4), has 
been known for a long time [7]. On the other hand, the temperature derivatives 
of the free energy, such as the internal energy or the specific heat, are highly 
nontrivial quantities which were obtained only in the stability region F < 2 of 
the 2D TCP by exploring the equivalent sine-Gordon model (for a short review, 
seeRef. [8]). 

The 2D OCP, which will be of interest in this paper, is formally related 
to the eigenvalue distribution of certain complex random matrices [9] and to 
the normalization problem of the Laughlin's wave functions in the fractional 
quantum Hall effect [10, 11]. There are indications from numerical simulations 
that, around F ~ 142, the fluid system undergoes a phase transition to a 2D 
Wigner crystal [12]. The existence of this transition has been put in doubt in 
a more recent paper [13]. As was already mentioned, by mapping onto free 
fermions the model is exactly solvable at the coupling F = 2, in the bulk [14] as 
well as in some inhomogeneous situations [5, 6]. The other exact information 
comes from the sum rules for truncated particle correlations valid at arbitrary 
F of the fluid regime. The usual zeroth- and second-moment conditions [1], 
having analogue in any dimension, are supplemented by the fourth-moment 
(compressibility) sum rule [15], available explicitly due to the knowledge of the 
exact equation of state, and the sixth-moment condition [16], related to universal 
finite-size properties of the Coulomb system [17]-[19]. At couphngs F = 27 (7 a 
positive integer), the partition function of the 2D OCP confined to a domain can 
be calculated exactly up to a relatively large finite number of particles N. For 7 
being an odd integer, the methods based on the expansion of even powers of the 
Vandermondc determinant into Schur functions [20, 21] turn out to be especially 
efficient. For 7 being an even integer, representations based on the permutation 
group [22, 23] are useful. At arbitrary integer 7, the 2D OCP is mappable onto 
a discrete ID fermionic field theory [24]. Within this fermionic representation, 
a symmetry of the model with respect to a complex transformation of particle 
coordinates has been shown to imply a functional relation for the two-body 
density. The functional relation is equivalent to an infinite sequence of sum rules 
relating the cocfBcicnts of the short-distance expansion of the two-body density. 
The generalization of the symmetry to multi-particle densities, possessing a 
specific invariant structure, was presented in Ref. [25]. 

The mathematical formulation of the 2D OCP looks at first sight simpler 
than the one of the 2D TCP. The mentioned integrability of the Coulomb gas 
therefore evokes the potential possibility of the integrability of the 2D OCP, and 
this is the main subject of the present paper. The integrability of the 2D jellium 
is investigated first at arbitrary coupling F via an equivalent 2D Euclidean-field 
theory, and then at special values F = 2*integer via the equivalent ID fermionic 
model introduced in Ref. [24] . As a by-product of the developed formalism, the 
gauge invariance of the 2D OCP is proven at the free-fermion point F = 2. 
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The paper is organized as follows. 

Section 2 is devoted to the 2D Euclidean-field representation of the 2D OCP. 

In Subsection 2.1., wc sketch shortly the phcnomcnological Dcbyc-Hiickcl calcu- 
lation of the free energy in order to have a test formula for functional methods. 
The 2D OCP is mapped onto a 2D Euclidean-field theory in Subsection 2.2. 
A comparison with the sinc-Gordon representation of the 2D TCP is made in 
Subsection 2.3. In Subsection 2.4., the "classical" integrability of the Euclidean- 
field representation of the 2D OCP is investigated by using a scheme proposed 
by Ghoshal and Zamolodchikov [26]. 

Section 3 is devoted to a further development of the discrete ID fermionic 
representation of the 2D OCP at couplings F = 2*integer [24]. At these cou- 
plings, the partition function of the plasma is shown to admit a representation 
in terms of a linear set of equations. 

Section 4 deals with gauge invariance of the bulk 2D OCP at coupling F = 2 
which has been proven previously by more standard methods in Ref. [27]. 
The alternative proof of gauge invariance presented here is related to the exact 
inversion of a class of infinite-dimensional matrices, which is of mathematical 
interest. 

A brief recapitulation is given in Section 5. 



2 2D FIELD REPRESENTATION 

2.1 Debye-Hiickel Calculation 

In the mean-field approximation, the effective electric potential tp at distance r 
of charge q, placed at the origin and surrounded by mobile g-charges plus the 
neutralizing background, is given by the 2D Poisson equation 

AV'(r) = -27rg {(5(r) + n [e-^'^^W _ ij | (2.1) 

The mean-field Boltzmann factor can be linearized for high temperatures. Eq. 
(2.1) then transforms to 

(A - K^) V(r) = -27rg(5(r) (2.2) 

where K is the inverse Dcbye length defined by = 27rrn. The solution of (2.2) 
reads V'(r) = qKo{Kr), where Kq is a modified Bessel function. The excess (i.e., 
over ideal) free energy per particle, /ex, is expressible in terms of the excess 
potential energy per particle 

Wex(/3) = |lim[V'(r)+9lnr] (2.3) 

z 1 — »o 



as follows 



/3/ex= / d/?'nex(/?') (2.4) 
Jo 
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Using the short-distance expansion of Kq [28], 

Ko{x) = -ln{x/2)-C + 0{x'^\nx) (2.5) 
C is the Euler's constant, we arrive at the expression 

/3/e.~-^ln(^^)+^(l-2C) (2.6) 

valid in the small coupling limit F — > 0. In what follows, this will be a test 
formula for functional methods. 

An analogous procedure can be applied to the F — > limit of the 2D TCP 
of ±q charged particles. The excess free energy per particle is again obtained 
in the form (2.6). 

2.2 Field-Theoretical Representation 

The 2D Coulomb potential (1.2) is singular at r = 0. This causes mathemat- 
ical difficiiltics when representing interacting Coulomb systems as equivalent 
field theories. To avoid this problem, we will consider the Coulomb potential 
regularized smoothly at short distances: 

Vreg(r) = - Inr - iiTo , e > (2.7) 

In 3D, the analogous regularization has been used in Ref. [29]. Since the Bcssel 
function Ko{x) decays to zero exponentially (r) has the large- 

r asymptotic of the pure Coulomb potential. On the other hand, using the 
short-distance expansion (2.5) in (2.7), the self-energy is finite 

t;reg(0) =C-ln2- ilne^ (2.8) 

It is easy to verify that the regularized Coulomb potential satisfies the diff'eren- 
tial equation 

(A - e^A^) vresir) = -27r5(r) (2.9) 

which is the counterpart of the 2D Poisson equation (1.1). 

We are interested in the bulk thermodynamic properties of the OCP, defined 
in the infinite region K = B? with the volume |A| oo. For N mobile particles 
at positions in we introduce the microscopic density of particles n 

and of the total charge p as follows 

N N 

n(r) = ^(5(r- r^), p{r) = q^5{r - Vj) - qrib (2.10) 
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Here, —qrib is the fixed (i.e., A^-independent) charge density of the background. 
The total interaction energy of the particle-background system is expressible as 



EN{{rj}) 



d r d r p{r)vi. 



|r-r'|)p(r')-iiVgVeg(0) (2.11) 



We will work in the grand canonical ensemble with the fixed background [30, 31] 

and position-dependent fugacity z(r) of particles. The grand partition function 
S at inverse temperature /3 is defined as the sum over all iV-particle states 



OO ^ 

y - 



JV=0 



N 



rjZ 



(r,)] 



-/3Bjv({r,}) 



(2.12) 



The multi-particle densities can be obtained as the functional derivatives of the 
generator S with respect to ^(r); after the functional derivatives are done, the 
homogeneous regime with the uniform fugacity z{r) = z is considered. At the 
one-particle level, the particle density is given by 



n = (n(r)) 

, A SE 



'ESz{r) 



(2.13) 



uniform 



At the two-particle level, one introduces the two-body density 



722 (r,r') 



(n(r)n(r')) - n5{r - r') 



= z{r)z{r') — 



S 5z{r)5z{r') 



(2.14) 



uniform 



The grand partition function (2.12) can be expressed in terms of a 2D 
Euclidean-field theory. We start by the standard procedure (see, e.g., Ref. 
[32]) and substitute the representation (2.11) of Ej^ in the Boltzmann factor 
exp(— /?i?jv)- The self-energy term rcnormalizes the fugacity, z(r) z{r) = 
2(r)exp[r?;reg(0)/2]. According to relation (2.9), -(A - e^A'^)/{2TT) is the in- 
verse operator of Vj-^g. The bilinear term in exp(— /JE'jv) can thus be linearized 
by applying the Hubbard-Stratonovich transformation 



exp 



j dV I dV'p(r)z;reg(|r-r'|)p(r') 

/ D^exp {/ d^r [i0(A - e^A^)^ + i,/2^ct>p] } 
/ V(j) exp [/ dV i0(A - £2 A2)0] 



(2.15) 



Here, ^(r) is a real scalar field with all derivatives vanishing at infinity and / Vcj) 
denotes the functional integration over this field. The terms (pAcp and (f>A'^(j) can 
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be turned into — |V0p and (At/))^, respectively, after performing integrations by 
parts with vanishing contributions from infinity. Inserting p from (2.10), particle 
coordinates in (2.12) become decoupled from each other and one can sum over 
N , with the result 



\\V4>? + \e\^<t>f - ^(r)e^^* + l^f2^n^,4> 



where 



D 



V(f)cxp <- d 



i|V0p + le2(A^)2 



(2.16) 

(2.17) 



is the normalization constant. In the homogeneous regime z{r) = z, the uniform 
shift in (p 

\n{z/nb) 



^(r) ^ ^(r) + 



iV27rr 



factorizes out the z-dependcncc of 



= exp 



|A|ln - +|A|nb 
rib 



(2.18) 



(2.19) 



X exp |- J dV i|V./>p + ie2(A0)2 + (-e'^^ + 1 + iV2nr<P^ | 

The particle density n is yielded by the homogeneous analogue of Eq. (2.13) as 
follows 

This means that from the grand partition sum (2.12) only the term with the 
strict system neutrality survives, in the spirit of Ref. [30]. The density- fugacity 

relationship is trivial, namely the density does not depend on the fugacity. This 
enables us to pass to the canonical ensemble via the Legendre transformation 



-l3F{n) 



In S - TV In z 



(2.21) 



where F is the Helmholtz free energy and N = n|A|. The excess free energy, 
related to F as follows — /3Fex = —/3F + Nlnn — N, then reads 



-/3Fex(n) = ^('c-ln2-ilne2 



Ini? 



(2.22) 



where we have substituted the explicit form of the self-energy (2.8) and grouped 
the field part into the quantity R defined by 



R 



/D<^exp { - /dV [i|V<^|2 + ie2(A<^)2 + „ (^_eiV2^0 + i + 1^2^^)] | 



/ D<^ exp { - / dV [i I 2 + ie2 ( a<^)2] } 



(2.23) 



For the pure Coulomb interaction (e = 0), the field representation of the free 
energy similar to the one described by Eqs. (2.22) and (2.23) was established 
directly in the canonical format by Brilliantov [33] . The problem of the divergent 
self-energy was incorrectly ignored there, although this one enters into the final 
formulae. In what follows, we aim at deriving the small-F expansion of Ini? in 
(2.22) in order to show that, in the limit e — > 0, the divergent self-energy term 
(X Ine^ is canceled and the Debye-Hiickel result (2.6) is reproduced correctly. 
For small F, we expand the exponential 

giv^^ _ ^ i^/2^(p - i(27rF)(/.2 (2.24) 

R then becomes equal to 

_ /P</>exp {- /dV [l|V0p + ^e^Acbr + 1^^'] } 

/C0exp{-/d2r [i|V<^P + ie2(A<^)2]} > 

The Gaussian functional integrals can be diagonalized in the Fourier k-space, 
with the result 

(2.26) 

The integration over k can be carried out explicitly and one finds 

\M 

Stt 
where 

1 ± VI - 4e2/t2 
t± = (2.28) 

In the e ^ limit, 

1 

e 

Consequently, 



lni?= (t+lnt+ + t-lnt- + \lne^ ) (2.27) 



t+ = ^ - + 0{e^), t-=K'^ + 0{e^) (2.29) 



jVF 

lni?~ — (Ine^ +lnK2 - 1) as e ^ (2.30) 

Inserting this into (2.22), the singular Ine^ term disappears and one recovers 
the Debye-Hiickel result (2.6). 

The two-body density can be obtained from the field representation (2.16) 
of S[2] using formula (2.14). The uniform shift in the (p-tie\d, relation (2.18), 
then leads to 

n2(r,r') (e'^^'!^('-)o'^^<^('-')) 



i(r)n(r') (^QiV2^4>{r)^e'-^^'^(''''>) 



(2.31) 
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where the averages (• • •) arc taken with the field action 



S[<P] 



l|V0p + le2(A.^)2+n(-e- 



+ l + i^ 



(2.32) 



Note that because the self-energy does not enter explicitly into (2.31), one can 
put e = in the action (2.32) (the consequent singularities in the numerator and 
the denominator must be precisely canceled with one another), and consider 



S[(j)]= j d^r i|V(/)p+n(-e'^^'^ + l+iV2rf^</)) 



(2.33) 



In the Dcbye-Hiickel F ^ limit, the expansion of the exponential according 
to Eq. (2.24) transforms the action (2.33) to 



'DH 



(2.34) 



Since (0) = with this action, the Wick's theorem for Gaussian integrals implies 



^giv^[^(r)+^(r')]^ = exp{-7rr([<^(r) + </)(r')]')} 



Consequently, 



^^^%^=exp{-27rr(0(r)0(r'))} 



(2.35) 
(2.36) 

(2.37) 



For the quadratic action (2.34), the correlator {(j){r)(j){r')) is equal to the inverse 
matrix element of the operator —V^ + k^, 



{^{r)^{r')) = -Ko{^\r-r'\) 



Thence 



712 (r,r') 



= e-r^o(''l'-'-'l) ~ 1 - TKo{K\r - r'\ 



(2.38) 



(2.39) 



in agreement with the standard Debye-Hiickel calculation (see, e.g., Ref. [16]). 



2.3 Comparison with the 2D TCP 

It is well known that the 2D TCP of ±q charges is equivalent to the sine-Gordon 
theory [32]. Namely, the grand partition function S is expressible as 



E{z) = 



J V<pex-p{-S{z)) 
JV(l>exp{-S{0)) 



(2.40) 
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where S is the 2D Euchdean sine- Gordon action 



S{z) = y" dV ^|V<;/)|2-2^cos(\/2^<;/)) 



(2.41) 



As before, z is the particle fugacity renormahzed by the (divergent) self-energy 
term, z = z cxp[/3u(0)/2]. In contrast to the previous case of the OOP, the 
self-energy v{0) cannot be eliminated from the functional integration through a 
uniform shift in 0. To give a precise meaning to z one has to fix the normalization 
of the coupled cos-field. In the Coulomb gas, the behavior of the two-body 
density for oppositely charged particles is dominated at short distance by the 
Boltzmann factor of the Coulomb potential. 



n+_(r,r') 



~ z'^lr ■ 



■r'|-^ aslr-r'l 







(2.42) 



which is the crucial supplement of the mapping (2.40) and (2.41). Under this 
conformal short-distance normalization, the divergent self-energy factor disap- 
pears from statistical relations established within the sine-Gordon formulation. 

One can document the above scheme in the Debye-Hiickel limit F — > 0, when 
cos(y2rf'0) ~ 1 - (27rr)(/.V2 and 



-\V(I>\'^ + 2'KTz(t>'^ 



- 2z\A\ 



(2.43) 



The total particle density and the two-body density of the opposite charges are 
given by 



n = 2z{e 



i\/27rr0(r) 



2 / i v^0(r) „-iy2^</.(r' ) \ 



z'{e 



(2.44) 
(2.45) 



respectively: the averages are taken with the action S'dh- The Wick's theorem 
for Gaussian integrals gives 



^giV2rf'[<^(r)-</.(r')]^ 

Simultaneously, 



exp{-7rr([</.(r)-<^(r')]2)} 



(.^(r).^(r')) = ^Ko (Vi^ir - r'|) 



with the short-distance asymptotic 



(</.(r)0(r')) 



1 

2^ 



In (V^|r-r'|) -C 



as r — r 



(2.46) 
(2.47) 

(2.48) 

(2.49) 
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Combining Eqs. (2.44)-(2.49) with the short-distance normahzation (2.42), one 
arrives at the expHcit density-fugacity relation 



— = {iiLY^^ /re 



2. \2J ^^H^J ^'-''^ 

valid in the leading F-order. This relation does not contain the self-energy v{0); 
the latter was coupled only to the power. Having the explicit n — z relation 
it is straightforward to pass to the canonical format. The Debye-Hiickel relation 
(2.6) is reproduced again. 

2.4 Classical Non-Integrability 

The field action of the 2D OCP (2.33) belongs to a more general class of actions 
possessing the local form 



S[4>] = I dV 



(2.51) 



where the factor 4 in the potential term appears for notation convenience. The 
exact solvability of a theory depends on the particular form of the potential V, 
which is in our case 

y(<^) = ^ (-e'^4> + 1 + (2.52) 

It is well known [34] that the 2D Euclidean field theory (2.51), defined in the 
space of points r = {x,y), is the imaginary- time {y = it) continuation of the 
equivalent real-time 1-1-1 dimensional quantum field theory with the action 



/OO POO 
dt / da; 
-CSO J — OO 



(2.53) 



This action is dominated by fields satisfying the "classical" equation of motion 
SS = 0. In terms of light-cone coordinates defined by d± = {dt ± dx)/2, the 
equation of motion reads 

d+d-(l) = -V'{(p) (2.54) 

In the case of our potential (2.52), rescaling appropriately the ^field, this equa- 
tion is nothing but the 1 + 1 analogue of the usual mean- field Poisson-Boltzmann 
equation for the OCP. In general, the integrability of a field theory is associ- 
ated with existence of an infinite sequence of conserved quantities (integrals of 
motion or "charges"). In what follows, we use a general scheme [26] to find out 
whether or not there exists an infinite sequence of conserved quantities for our 
field theory with the potential V given by (2.52), in the classical limit, i.e., when 
the field (f> is governed by the equation of motion (2.54). The scheme represents 
a unique way of determinig integrability properties of the given field theory. 
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Existence of a conserved charge is associated with the appearance of a pair 
of "conjugate" local field densities {T,0), with zero boundary conditions at 
a; — > ±00, such that 

d-T = d+e (2.55) 

In terms of x and t variables this is equivalent to dt (T — 9) = dx {T + 0). The 
integration over x results in 



d_ 
dt 



/OO 
da;(r - 61) 
-CSO 



iT + 9)r=0 (2.56) 



and thence the charge / dx{T — d) is conserved. We look for (T, 9) as polynomial 
functions in derivatives 9+0, d'^(j), etc. The notation T, {9s) will be used for T (9) 
with just s 5+-derivatives. Because of the specific form of the equation of motion 
(2.54), only T^+i and 9s-i can create the conjugate couple, d-Tg+i = d+9s-i- 
The conserved charge is then 



Qe = j Ax{Ts+i-9s-i) s=l,2,... (2.57) 

Note that the total d+ derivatives can be dropped from T because the consequent 
difference d+ — d- = dx produces vanishing boundary contributions to the 
conserved charge. Similarly, let Tg-i and 9s+i be the conjugate polynomials of 
given orders in 9_ -derivatives such that d-Tg-i = d+6s+i- Then, the charge 

Qs^ j Ax{9s+i-fs-i) s=l,2,... (2.58) 

is conserved. As before, the total 9_ derivatives can be dropped from 0. 
At s = 1, writing T2 = {d^<j))^ one has 

d-T2 = 2{d+<j>){d+d-(t)) = -2{d+(^)V'{(^) = d+ [-2V{(t))] (2.59) 

so that ^0 = — and 

Qi = j dx [{d+(t>f + 2V] (2.60) 

Analogously, h = {d-cpf, Tq = -2V and 

Qi= J dx[{d-,^f + 2V] (2.61) 

Qi + Qi and Qi — Qi are energy and momentum, respectively, and these two 
quantities are always conserved for any potential V. There is no solution for 



12 



conjugate polynomials producing conserved charges at s being an even integer. 
For s = 3, there exist conjugate polynomials 

T4 = (^^^\d+<i^r+{di<j>r (2.62) 

02 = -{d+(t>fv"{(l>) (2.63) 

and the corresponding {64, T2), provided that the potential satisfies the differ- 
ential equation V'" = b^V' . This equation is fulfilled either for the trivial free 
field theory {b = 0) or for the potential 

V{<f)) = Ae'"'' + Be-'"'', 6^0 (2.64) 

When the constants A and B are nonzero and equal to one another, one recog- 
nizes the sinh-Gordon (6 real) or sine-Gordon {b imaginary) models. At s = 5, 
there exist conjugate polynomials and the corresponding conserved charges if 
the potential is either of the previous form (2.64) or of the form 

y(^) = Ae'"'' + Be-(*/2)"^, b^O (2.65) 

This integrable field theory is known as the BuUough-Dodd model [35] and for 
imaginary b it corresponds to the 1 : 2 charge-asymmetric Coulomb gas [36]. 
The models with potentials (2.64) and (2.65) arc the only two one-component- 
field members of the integrable afiine Toda field theories, based on the Dynkin- 
diagram classification of simple Lie groups. Programming the whole scheme 
in the symbolic language Reduce, we were able to proceed up to the relatively 
high s = 15 order. Except for the repeated appearance of the two potentials 
(2.64) and (2.65), we did not find any other solution for the potential leading to 
conserved charges. We do not anticipate a sudden appearance of an additional 
potential producing conserved charges for s > 15. If it is so the 2D OOP, 
characterized by the field potential (2.52), is not classically integrable. 

The conjectured classical non-integrability does not exclude the complete 
"quantum" (all realizations of the (/)-field are considered) integrability of the 
model. At specific values of the coupling constant F, quantum fluctuations of 
the field around its classical saddle-point value can make the plasma integrable, 
as it is at the free-fermion point F = 2. 



3 LINEAR FERMIONIC REPRESENTATION 

We now consider the 2D OCP, confined to a domain A, directly in the canoni- 
cal ensemble. N pointlike g-charges are embedded in a spatially homogeneous 
background of charge density pi, = —qrib. If the system-neutrality condition is 
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imposed it holds Ub = N/\A\. The background produces the one-body electric 
potential Vb{r) = pb /^d^r't;(|r — r'|) which satisfies the Poisson equation 

Avb{r) = 2TTqnb, r e A (3.1) 

Since, written in the complex {z, z)-coordinates, 

A = Ad,d, (3.2) 

for a circularly symmetric background one has 

irqnb _ 

fcirc = const H — zz (3.3) 

The deformation of the circular boundary dA or the presence of some charge 
densities outside of A generates an additional gauge potential Ugauge such that 
Atigauge(-z, ^) = for 2; e A. With regard to (3.2), the general solution of this 
equation reads 

oo 

%auge(2;, z)=Ao + Y^ (A.z' + B,z') (3.4) 

s=l 

Since Coulomb potentials are real, physical situations correspond to the choice 
As = Bg for all s. The particular case of a quadrupolar potential Vgauge = 
A{z'^ + z^) results from the deformation of the disk into an ellipse [11, 37]. 

The potential energy of N particles at positions {zi € A} plus the back- 
ground is 

E = Eq + q^Vb{zi, Zi) -q^^\n\zi- zj \ (3.5) 

i i<j 

The background-background interaction constant Eq does not influence the par- 
ticle densities and so it can be omitted. The partition function at inverse tem- 
perature /3 reads 



1 f ^ 

■^Af = / n ['i'^Z,w(Zz, Z^)] Yl I 

■ •'^ i=l i<j 



N 

\zi-zjf (3.6) 

i<j 



where F = Pq^ and the one-body Boltzmann factor w{z, z) = eyiY)[—(}qVb{z, z)]. 
The multi-particle densities can be obtained in the standard way [see relations 
(2.13) and (2.14)]: 

n{z,z) = w{z,z)^ Jf^ (3.7) 
Z]S! bw\z^z) 

ni{zx-,zx\zi,Z2) = w{zi,zi)w{z2,Z2)-^-^, — 1 f ^, ^ (3.8) 

Zn dw{zi,Zi)w{Z2,Z2) 

etc. Wc will study the special case of the plasma with a "soft wall" [38] when 
particles are confined by the background itself. In particular, for a fixed value 
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of the particle number N one makes the radius of the homogeneous circular 
background infinite, |A| — > oo. The particles will gather in a circular region of 
area N/rib. In the limit N ^ oo, the soft- wall edge of this region will have the 
structure different from the one of the usual hard- wall problem, but the particle 
densities deep in the interior (|^;| finite) will not be touched by the soft wall 
and will attain the bulk values. Also the free energy will be modified only by 
a surface term. Within the soft-wall |A| ^ oo formulation of the problem, we 
will consider two cases: 

(i) the circularly symmetric homogeneous background, Eq. (3.3), with 

1 / r _\ 
w{z, z) = — cxp ( — —TTUhZz 1 (3.9) 

(ii) the general homogeneous background with the gauge component (3.4), 
given by 

p oo 

■u;(z,z)=cxp ——Trni,zz — ''^{asz'' + bsZ'^) (3.10) 

s=l 

where Ug = PqAs and bg = PqBg. As was already mentioned above, real physical 
situations require Ug = bg for all s and the values of gauge parameters must be 
such that the multiple integral determining the partition function (3.6) does not 
diverge. 

For the coupling F = 27 (7 an integer), it has been shown in Ref. [24] 
that the partition function (3.6) can be expressed as the integral over two sots 
of Grassmann variables {^i°'\4'i"^} each with 7 components (a = 1,...,7), 
defined on a discrete chain of A'' sites i = 0,1, . . . ,N — 1 and satisfying the 
ordinary anticommuting algebra [39], as follows: 

ZAr(7) = J VtpV£, e^^^'^'^ (3.11) 

7(JV-1) 

5(C,^) = E ^i^ij'^o (3-12) 

Here, VijjVS, = riilo^ dTji'-'*'^ . . . d^'-'^^d^^'''^ . . . d^|^^ and S involves pair interac- 
tions of "composite" operators 

E E (3.13) 

«!,■■■, «7 = Jl , ■ ■ ■ , J7 =0 

(»l + ---+»7 = ») Ul + ---+3^=j) 

i.e., the products of all 7 anticommuting-fiekl components, belonging to either 
^- or ip-set, with the fixed sum of site indices. The interaction strength is given 

by 

Wij= [ d^zw{z,z)z'z^; i,j = 0,l,...,^{N-l) (3.14) 
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Using the notation (•••)= / VipVS^e^ ■ ■ ■ jZ^i^l) for an averaging over the anti- 
commuting variables, the particle density (3.7) and the two-body density (3.8) 
are expressible in the fermionic format as follows 

7(JV-1) 

n(z,z) = w{z,z) ^ (5i*j)z*z^ (3.15) 

n2{zi,Zi\Z2,Z2) = w{zi,Zi)w{Z2,Z2) 

X (^i^^j^iJ^j^z^z^'z^z^ (3.16) 

n.jii«2 ,32=0 

respectively. 

The exact solvability of the 2D OCP at F = 2 (7 = 1) is due to the bilinear 
form oi S = J^ljlo^^ ^iWijtpj. Thus, 

Z;v(7 = 1) = ^et{wij)\lrj^ (3.17) 

The two-correlators determining the particle density (3.15) are equal to the 
inverse elements of the N x N w-matrix (3.14), 

i^ii'j) = wj' (3.18) 
The Wick's theorem applied to the four-correlators in (3.16) implies 

= ^7X^j2i ~ ^j2i'^jX (3-19) 

In the case of a circularly symmetric plasma with w{r) = w{r) and A = {r < 
R} [like the one of interest, defined by Eq. (3.9) and i? — > 00], the interaction 
matrix w with elements (3.14) becomes diagonal: 

R 

2i+l 



WiSij, Wi = 2-K / d™(r)r'''+^ (3.20) 



The "diagonalized" form of the partition function 

„ 7(A'-l) 

Zjv(7) = / 'D^'Di n exp(Si«;i*i) (3.21) 

i=0 

implies that only correlators {^i^'^ j^'^i^^ ■ ■ ■) with 11+12 + ■■ ■ — ji+ 32 + ■ ■ ■ 
will be nonzero. The dependence of .^jv(7) on the set of moments {wi):1^~^'' 
is the crucial problem whose solution would mean the complete exact solution 
(free energy, correlation functions) of the bulk 2D OCP at the given coupling 



16 



r = 27. Let us write down explicitly few examples of this dependence. At 
7 = 1, we have the simple result 

Zn{1)=woWi---wn-i (3.22) 

At 7 = 2, using the anticommuting integral rules one finds from (3.21) for small 
particle numbers N = 2,3 that 

^2(2) = woW2 + 2wl (3.23) 
^3(2) = W0W2W4, + 2wowl + 2wl'W4 + 4'Wi'W2W3 + 6W2 (3.24) 

etc. At 7 = 3 one has 

^2(3) = woW3 + 3'^wiW2 (3.25) 
■^3(3) = W0W3WQ + 3^woW4W5 + 3^wiW2We 

+6'^wiW3W5 + 15^u'2W3W4 (3.26) 

etc. There exists one model exactly solvable for every 7 and N, namely the 
2D OCP constrained to a circle. In that case w{r) = S{r — l)/(27r) and, conse- 
quently, = 1 for alH = 0, 1, . . . , 7(A'' — 1). It was proved in various ways [40] 
that 

hNV. 

Relations (3.22)-(3.26) pass this test. 

We now aim at analyzing the structure of a general summand in Zn{j). 
It follows directly from the fermionic representation (3.21) that each term is 
composed of just N w's, Wi-^Wi^ ■ ■ -Wij^. The transformation A^^^"-* for 

all a = 0, 1, . . . , 7 indices and sites i = 0,1, . . . , N — 1 implies S A'S. As 
a consequence, the subscripts of the general term Wi^Wi^ ■ ■ -Wij^ must satisfy 

the relation ^ i = 7(0 + IH V N - I) = -fN{N - l)/2. It is necessary to 

distinguish between 7 being an odd or even integer. 

For 7 an odd integer, the composite S and ^' operators are products of an 
odd number of anticommuting variables and so they themselves satisfy the usual 
anticommutation rules 

{Ei,E,} = = {S„vp,} = (3.28) 

for all i,j = 0,1,..., j{N - 1). In particular it holds = = q. The 
expansion of each exponential in (3.21) is thus 

exp{EiWi^i) = 1 + EiWi-^i (3.29) 

and the partition function is represented as follows 

Zn{i) = C'il,i2,...,iiv^ii^^2 • • •W^^iv (3-30) 

il<i2<-<ijv=0 
[^i=-f«(iV-l)/2] 
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Cil,i2v>»N 



(3.31) 



We see that for odd 7 a given Wi can occur in a summand at most once. 

When 7 is an even integer, the composite S and ^' operators are products 
of an even number of anticommuting variables and so they commute with each 
other 

[E„E,] = [%,^j] = [E,,^^]=0 (3.32) 

for all i,j = 0, 1, . . .,'y{N — 1). The expansion of an exponential is now a bit 
more complicated, 



EfX''^^' jii^im^iV for i = 0, . . . , 2(7V - 1) 

[ ^^2f-i-PV7] fori = 2(Ar-l),...,7(Ar-l) 



(3.33) 

Next terms vanish because they contain second or higher power of at least one 
anticommuting ^ or tp variable. Let us introduce the indices {ai}]l^~^^ with 
the following value ranges: 



0,1,..., 



0,l,...,2Ar- 1 



fori = 0,...,2(iV-l) 



fori = i(Ar-l),...,7(7V-l) 



(3.34) 



The partition function for 7 an even integer is then expressible as 

an ai ^-f{N-l) 



7(JV-1) 



where 



oO.'-',<»-,(JV-l) 



ao,ai,...,a,(«-i) ^olai ! • • • a^(iV-l) ! 



TIC ■^"OTT"! 
'^'i "0 "1 



„Q^(JV-1) 



(3.35) 



(3.36) 



The underlying fermionic representations, relations (3.30), (3.31) for 7 odd 
and relations (3.34)-(3.36) for 7 even, contain the unknown coefficients which 
can be formally written in both cases as 



C'n,...,ijv — / ^n^i2 



The "basic sector" of indices is 



ll < l2 < ■ ■ 
il < 12 < ■ ■ 



< iN, for 7 odd 

< In, for 7 even 



(3.37) 



(3.38) 
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The C-coefRcient with an arbitrary sequence of indices can be expressed in terms 
of the basic one with indices ordered according to (3.38) through a successive 
exchange of nearest-neighbor couples of composite operators, 

= (~l)'''Cii,...,ij_|.i,ij,...,ijv (3.39) 

It is evident that 

Co,7,...,7(JV-i) = 1 (3.40) 

One can cvahiatc the C-cocfRcicnts directly from their definition (3.37) by 
using the rules of integration over anticommuting variables. However, there 
exists another simpler way how to determine these coefficients. Let Pj-i{k, k') 
be a general polynomial of the (7 — 1) order in the k, k' variables, and consider 
this polynomial in the following combination with the composite operators 

7(JV-1) 

J2 Pj-iik,k')EkEkr, K = 0,l,...,2^{N-l) (3.41) 

k,k'=0 
{k + k' = K) 

Representing the S-operators in terms of the anticommuting ^-variables [see 
relation (3.13)], Eq. (3.41) can be rewritten as 

J2 P^-i{ki+- ■ ■+k^, k[+- ■ ■+k'^)i'^ ■ ■ ■ Ci^^^g) • • •Ci7)<5K,fe,+...^+fe;+...^ 

*=! '=7 _Q 

(3.42) 

Since Pj-i is the polynomial of the (7 — 1) order in its arguments, at least 
one couple of the k, k' factors with the same subscript, say the couple ki and 
k[ associated with does not occur in the given term of the expansion of 
P^-i{ki -\- ■ ■ ■k^,k'iA k'^). Thus, since 

JV-l 

Et(i)t(i)A 

fci,fei=0 

= E (^i'^^? + ^il^^if ) 5KM+-k.+K+-k', = (3.43) 

the polynomial combination with composite operators (3.41) vanish. Passing 
from F-y_i(fc, k') to Q^-i{k + k',k — k') one has 

7(Ar~l) 

J2 (fc - fc')"^fcSfc' = forn = 0,l,...,7-l (3.44) 

k,k' = 
(k + k' = K) 
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We remind that the S-operators anticommute for 7 odd and commute for 7 
even. This is why for 7 odd only odd powers of (k — k') provide a nontrivial 
information and for 7 even only even powers of (fc — k') arc informative. In view 
of the representation (3.37), we finally arrive at a homogeneous set of linear 
equations for the coefBcients C: 
if 7 is odd, 

J2 (fc-fc'rC'n,...,iiv-2,fc,fc' =0; n= l,3,...,7-2 (3.45) 

k,k'=0 
{k-\-k' = K) 

if 7 is even, 

7(iV-l) 

J2 {k-kTCi„...,^^_,,k,k'=0; n = 0,2,...,7-2 (3.46) 

k,k' = 
(k + k' = K) 

In both cases K runs over 0, 1, ... , 27(iV — 1). 

In the first nontrivial case of 7 = 2, according to Eq. (3.39) the interchange 
of subscripts does not change the sign of the C-coefficients. The normalization 
is C'o,2,...,2(Ar-i) = 1- Formula (3.46) with n = applies. For N = 2, the 
normalization Co, 2 = 1 is supplemented by the only equation 

= Co,2 + + C2,o (3.47) 

so that Ci,i = —2. Inserting the two coefficients into the representation (3.35) 
one recovers the exact relation (3.23). For = 3, the normalization Co, 2,4 = 1 
is supplemented by the (overcomplete) set of five equations 

= Co, 2, 4 + Co, 3, 3 + Co, 4, 2 

= Ci,i,4 + Ci,2,3 + Ci,3,2 + Ci,4,i 

= C2,0,4 + C'2,1,3 + C2,2,2 + C'2,3,1 + C2,4,0 (3.48) 

= C3,o,3 + (^3,1,2 + C'3,2,1 + C3,3,o 

= C4,o,2 + ^4,1,1 + C4,2,0 

which implies Co,3,3 = Ci,i,4 = —2, Ci,2,3 = 2 and C2,2,2 = —6. Inserting these 
coefficients into (3.35) one recovers the relation (3.24). We suggest that the sets 
of linear equations, (3.45) for 7 odd and (3.46) for 7 even, constitute, together 
with the relation (3.39) and the normalization (3.40), complete (more precisely, 
overcomplete) sets to be solved for the C-coefiicicnts. We have checked this 
conjecture for 7 = 2 (F = 4) up to iV = 13 particles and for 7 = 3 (F = 6) up 
to A'' = 9 particles. 

In comparison with the technique of anticommuting variables, the algorithms 
for solving a set of linear equations are much faster. The presented scheme is 
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therefore very convenient for exact computer calculations of the plasma with a 
finite number of particles. Prom the methods dealing with the finite number 
N of particles [20]- [23], the closest one [23] also provides the representations of 
the partition function (3.30) for 7 odd and (3.35) for 7 even. The C-coefhcients 
are expressed there as multiple integrals over the unit circle, and consequently 
as multiple summations over all possible permutations of N indices. Such algo- 
rithm is inferior to the present approach. 

At this stage, we were not able to find a simplification of the underlying sets 
of linear equations in the thermodynamic limit. In spite of this the fact that 
the crucial C-coefhcients are determined by linear equations may indicate the 
exact solvability of the 2D OCP at F = 2*integer. 

4 GAUGE INVARIANCE 

In the case of classical statistical systems with short-ranged interactions among 
constituents, the thermodynamic limit of an intensive quantity docs not de- 
pend in general on the shape of the system and on the conditions at the 
boundary given by the surrounding medium. This may be no longer true for 
macroscopic systems with long-ranged interactions. A typical example is the 
shape-dependence of the dielectric susceptibility tensor for Coulomb conductors 
[41] caused by the long-range decay of charge correlations along the boundary 
[42, 43]. Gauge invariance of Coulomb fluids is another, in a certain sense oppo- 
site, phenomenon related to the long-range nature of particle interactions. Let 
us consider a macroscopic Coulomb conductor in an arbitrary shaped domain 
A, with perhaps some charge distribution on the boundary OA. The effect of 
the boundary is then reflected in the bulk through a long-ranged one-body po- 
tential whose Laplacian is zero inside the domain A. The assumption of gauge 
invariance states that this one-body "gauge" potential is perfectly screened by 
the Coulomb system at macroscopic distances from the boundary, and so it does 
not affect the averaged particle distributions in the bulk interior. 

The above developed fermionic formalism will be now used to prove gauge 
invariance of macroscopic Coulomb fluids at the special value of the coupling F = 
2. The original proof of gauge invariance (up to the gauge potential being the 
polynomial of degree 2 in spatial x and y components) at F = 2 was presented in 
Ref. [27] for a class of possibly inhomogeneous backgrounds. The present proof 
is an alternative one, applicable also to non-physical situations with a complex 
gauge potential. 

At F = 2 and in the units of Trrif, = 1, the one-body Boltzmann factor of the 
circularly symmetric background (3.9) defined in an infinite space reads 

w{z,z) = -e-'^, K = E? (4.1) 
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The N X N matrix (3.14) becomes diagonal, Wij = WiSij, with elements 

r°° 2 

Wi = 2 dre"'' r^'+^ = i!; i = 0,1, . . . , N - 1 (4.2) 
Jo 

Combining Eqs. (3.15) and (3.18), the latter written as {^ii^j) = Sij/wi, the 
particle density at distance r from the center is given by 



n{r) 
rib 



E 



(4.3) 



In the limit N ^ oo and for finite r, the bulk particle density is constant, 
n(r) = rif,, as was expected. 

Let us now consider the general case of a homogeneous background charac- 
terized by the one-body Boltzmann factor (3.10) with gauge degrees of freedom. 



iu{z, z) = — exp 

TT 



-zz-^{asz'' + bsz'') 



The interaction (3.14) is non-diagonal. 



fd^z 



z^z^ exp 



-zz - ^(0^2:** + bsZ') 



(4.4) 



(4.5) 



i,j = 0,l,...,N — l. The strength of gauge parameters {as,bs} must be such 
that integrals in (4.5) converge. For a real gauge potential with Ug — bs = ae"^, 
writing z = re'^ we have agz'^ + bgZ^ = 2ar^ cos{(j) + sip). If s > 3, for a 
fixed there always exist the (^-angles such that acos(0 + s(p) < 0, which 
implies divergent r ^ oo contributions to the integral in (4.5). Thence all gauge 
coefiicients {as,bs} with s > 3 must be equal to zero in physically acceptable 
situations. For s = 2, the convergence is ensured provided that 2a < 1. 

In view of Eqs. (3.15) and (3.18), the particle density at point {z, z) is given 

by 

n{z, z) 



N-l 

= -KW{Z, Z) ^ 
i,i=0 



Wjlz'z^ 



(4.6) 



In the thermodynamic limit N ^ oo, the bulk particle density is expected to 
be constant, n{z, z) = n^. This is true iff 



i,j=0 



exp 



^(a^z^ + 6s^") 



(4.7) 



This equation can be understood as the generating relation, with z and z as- 
sumed as independent variables, for the inverse elements of the infinite w-matrix 
with elements (4.5). 
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To prove that the quantities w^^ determined by Eq. (4.7) form indeed the 
inverse matrix of the one with elements (4.5) we have to show that it holds 

oo C30 

X] wikw-^j = 6ij, J2 '^"zk^kj = Sij (4.8) 

fe=0 fe=0 

for all i,j = 0, 1, . . .. Let us consider the auxiliary function 

oo 

Ki{t) = mkw^-t^ (4.9) 

i,/c=0 

Presuming the validity of the generating Eq. (4.7), a series of straightforward 
transformations yields 

2, 



Ki{t) = J2 -^z'z''w{z,z)w^lt^ 



j,k=0 



n 



I 



z exp 

d^z 



-zz + - ^ as{z^ - t^) 



(4.10) 



TT 



z + ty exp < -zz -tz-^tts [{z + ty - t^] 



Since {z + ty — t'' = z^ + stz" ^ + • • • + st* ^z, only terms t and —zz survive 
from {z + ty and the exponential, respectively. Consequently, 

Ki{t)=f foralH = 0,l,... (4.11) 

With regard to the definition (4.9) of Ki{t), this equation implies the first re- 
lation in Eq. (4.8). The second relation can be obtained in a similar way by 

showing that Li{t) = X^j,fc=o ^ife^^fcj^'' is equal to ? for all i = 0,1, The 

proof is accomplished. 

One can readily show that the proof of gauge invariance for the one-body 
density automatically ensures gauge invariance for the two-body density (3.16) 
with correlators (3.19). 

Strictly speaking, gauge invariance has meaning only for real gauge poten- 
tials; the short discussion after formula (4.5) tells us that only gauge potential 
being polynomial of degree 2 in z and z complex coordinates is allowed in physi- 
cal situations. On the other hand, the proof of the matrix inversion (4.8)-(4.11) 
requires only the finite values of matrix elements (4.5), without putting any 
further restriction on the {as,bs} gauge parameters which may therefore be 
unphysical. Typical unphysical examples leading to finite values of matrix el- 
ements arc bs = or a., = bg = i. It is interesting from a mathematical point 
of view that there exists a large family of infinite matrices, with elements Wij 
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{i, j = 0,1, . . .) defined by Eq. (4.5), which are explicitly invertible by using the 
closed- form generating formula (4.7). A detailed analysis of this mathematical 
peculiarity goes beyond the scope of the present paper. 

5 Conclusion 

In the canonical ensemble, the 2D jcUium is equivalent to the 2D Euclidean-field 
theory with the action (2.33). Here, the divergent Coulomb self-energy does not 
renormalize the model's parameters like it is in the sine-Gordon representation 
of the 2D Coulomb gas. In contrast to the sine-Gordon model, the quantum 
analogue of the present Euclidean theory is conjectured to be not integrable 
on the classical level (only such realizations of the ^field are considered which 
minimalizc the action) due to the lack of an infinite sequence of integrals of 
motion. The classical non-integrability does not exclude the complete quantum 
(all realizations of the (p-&eld are considered) integrability at specific discrete 
values of T, like it is at the frcc-fermion F = 2 point. This frcc-fermion coupling 
belongs to a family of couplings F = 27 (7 integer) which admit a ID fermionic 
representation of the partition function. The fermionic representations, relations 
(3.30) and (3.31) for 7 odd and relations (3.34)-(3.36) for 7 even, contain the 
unknown C-coefficients. These coefficients are determined by the homogeneous 
sets of linear equations, (3.45) for 7 odd and (3.46) for 7 even, supplemented 
by the exchange formula (3.39) and the normalization (3.40). This feature is a 
sign of integrability. The present analysis might be a challenge for specialists in 
the Field Theory. 

The proof of gauge invariance of the 2D OCP at F = 2 is related to the exact 
inversion of a class of infinite-dimensional matrices which elements are deter- 
mined by non-Gaussian integrals (4.5). This is interesting from a mathematical 
point of view. 
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